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ABSTRACT 

In  this  report,  techniques  generally  employed  in  the  analysis  of  intercept  guidance  problems  are 
reviewed.  From  the  governing  non-linear  equations  describing  such  problems,  two  basic  linear 
models  are  derived.  Traditionally,  these  linear  models  are  utilised  as  a  basis  for  preliminary 
intercept  engagement  studies.  Under  certain  input  conditions,  the  two  models  are  mathematically 
equivalent  and,  hence,  have  been  used  interchangeably  by  weapons  analysts  to  yield  appropriate 
design  and  performance  data  in  support  of  their  programs.  However,  for  a  specific  set  of  initial 
conditions,  which  includes  a  very  important  class  of  practical  problems  that  may  be  assessed  with 
the  use  of  these  models,  it  is  noted  herein  that  one  of  these  linear  models  produces  incorrect 
performance  data  when  compared  to  a  non-linear  simulation  of  the  engagement.  In  contrast,  the 
other  model  produces  consistent  results  with  those  generated  by  the  non-linear  simulation 
regardless  of  the  initial  conditions  considered.  To  remedy  this  discrepancy,  the  necessary 
mathematics  are  derived  to  bring  the  two  formulations  into  alignment  for  any  form  of  the  initial 
conditions  and  inputs  to  the  system.  Consequently,  this  leads  to  a  consistency  in  the 
corresponding  adjoint  models  which  are  constructed  from  these  linear  models,  thus  ensuring  the 
generation  of  correct  output  data  regardless  of  which  model  is  employed  by  the  analyst. 
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On  Alternative  Formulations  for 
Linearised  Miss  Distance  Analysis 


Executive  Summary 


In  this  report,  the  missile-target  engagement  problem  is  analysed.  As  part  of  the  analysis, 
the  non-linear  governing  equations  are  reviewed  for  motion  in  a  single  plane.  These  non¬ 
linear  equations  are  employed  for  two  purposes.  Firstly,  the  equations  form  the  basis  of  a 
non-linear  simulation  program  developed  in  Simulink.  This  Simulink  model  is  controlled 
and  executed  via  a  graphical  user  interface  specifically  developed  as  an  aid  for  the 
weapons  analyst  to  study  the  engagement  problem.  Secondly,  the  non-linear  equations  are 
utilised  as  a  basis  for  linearisation  and,  hence,  the  derivation  of  an  approximate  linear 
model  of  the  engagement.  Two  different  formulations  of  the  linear  model  are  derived. 
These  are  designated  as  Model  A  and  Model  B  in  the  report. 

Although  both  models  are  generally  used  interchangeably  in  the  literature  for  guided 
missile  homing  loop  analysis,  it  is  demonstrated  herein  that,  under  certain  input 
conditions,  care  needs  to  be  exercised  when  using  one  of  these  models.  Model  A,  for 
performance  analysis.  In  order  to  ensure  that  both  models  yield  the  same  performance 
data  for  all  input  conditions  considered,  a  correction  factor  is  derived.  This  correction 
factor  needs  to  be  included  in  the  form  of  an  added  initial  condition  on  one  of  the  states  in 
the  state  space  representation  of  Model  A.  Simulation  results  show  that  the  two  models 
are  in  agreement  when  this  correction  factor  is  applied.  Knowledge  of  this  fact  is 
important  to  ensure  that  analysts  generate  correct  performance  data  when  using  linear 
techniques  such  as  the  adjoint  method.  The  adjoint  model  is  constructed  from  a 
knowledge  of  the  forward  linear  model,  that  is.  Model  A  or  Model  B,  and  is  traditionally 
employed  by  analysts  as  part  of  the  solution  process. 

To  gain  further  insight  into  the  nature  of  the  missile-target  engagement  problem  and  the 
parameters  that  influence  performance,  also  included  in  this  report  is  an  analytical 
treatment  of  the  problem.  It  is  well  known  that,  when  the  missile  guidance  dynamics  is 
represented  by  a  first  order  lag  (single  time  constant  system),  then  the  linear  differential 
equations  describing  the  intercept  problem  are  readily  amenable  to  analytical  treatment. 
Consequently,  an  analytical  investigation  of  each  model  (A  and  B)  is  carried  out  for  the 
range  of  input  conditions  considered  herein.  The  resulting  closed  form  solutions  from  each 
model  are  compared  and  are  shown  to  be  mathematically  equivalent  for  all  input 
conditions  considered  provided  that  Model  A  has  the  proposed  correction  factor 
implemented.  For  simulation  purposes,  the  corresponding  adjoint  model  of  Model  A  and 
Model  B  are  constructed  and  then  implemented  in  Simulink.  As  with  the  analytical  results, 
outputs  from  the  two  models  are  shown  to  be  in  agreement.  Finally,  the  models  are 
extended  to  represent  more  realistic  guidance  system  dynamics  (fifth  order  system)  and 
simulation  results  are  generated  and  compared. 
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Finally,  a  special  relationship  linking  two  of  the  derived  miss  distance  formulas  is  noted 
and  explored  further.  This  relationship  highlights  a  connection  between  the  miss  distance 
due  to  an  initial  target  displacement  and  that  due  to  an  initial  heading  error  in  the  context 
of  the  linear  analysis.  Following  verification  using  linear  simulation,  a  formula  based  on 
this  relationship  is  derived  and  proposed  as  a  means  for  predicting  performance  data  of 
more  complex  linear  systems.  This  formula  may  also  be  used  in  connection  with  the  non¬ 
linear  simulation  model  for  generating  approximate  miss  distance  profiles  due  to  the 
effects  of  a  step  in  target  position  prior  to  intercept.  A  step  in  target  position  typically 
arises  in  problems  associated  with  the  seeker  resolution  of  the  target  in  a  multi-target 
scenario.  It  may  also  arise  in  the  case  of  a  single  target  scenario.  For  this  case,  the  missile 
may  be  on  a  collision  course  with  the  predicted  intercept  point  (PIP)  of  the  target. 
Flowever,  at  the  time  of  seeker  turn-on,  the  PIP  may  not  necessarily  coincide  with  the 
actual  target  position. 
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1.  Introduction 


The  equations  describing  the  dynamics  and  geometrical  interactions  between  an  interceptor 
and  its  target  are  generally  non-linear  and  complex  and  are  usually  beyond  the  scope  of 
analytical  investigations.  Therefore,  they  are  solved  approximately  using  computer 
simulation.  Simulation  provides  the  missile  analyst  with  a  basic  capability  for  assessing  the 
performance  and  behaviour  of  the  interceptor  under  varying  conditions  of  the  engagement 
and  can  provide  answers  to  questions  such  as,  "What  is  the  overall  effect  on  miss  distance  due 
to  a  sudden  shift  in  target  direction?"  However,  to  gain  critical  insight  into  the  nature  of  such 
interactions  and  to  identify  key  parameters  that  may  affect  the  performance  and  behaviour  of 
the  interceptor,  engineers  often  make  use  of  simplifying  assumptions  in  order  to  linearise  the 
governing  equations.  Once  in  linear  form,  the  engagement  equations  are  more  easily  tackled 
using  established  linear  techniques.  Two  such  techniques  are  the  adjoint  simulation  method 
and  the  covariance  analysis  method. 

In  this  report,  we  first  review  the  engagement  model  based  on  two  dimensional  (2D) 
non-linear  dynamics.  Following  this,  we  make  use  of  linearisation  to  reduce  the  equations  to  a 
form  amenable  to  linear  analysis.  Using  standard  block  diagram  algebra,  these  linear 
equations  are  then  used  to  derive  two  basic  models,  referred  to  as  Model  A  and  Model  B, 
which  may  be  used  as  the  basis  for  the  linear  analysis  of  planar  missile/ tar  get  engagements. 
In  the  literature,  these  models  are  often  used  interchangeably  for  preliminary  analysis  of  the 
performance  of  guided  missile  systems  [1-6].  However,  as  is  demonstrated  herein,  caution 
must  be  used  when  adopting  one  of  these  models  under  certain  initial  conditions  of  the  state 
variables.  It  is  shown  that  the  mathematical  equivalence  of  the  two  block  diagram  topologies 
representing  the  models  is  dependent  on  the  nature  of  the  initial  conditions  imposed  on  the 
states  of  the  system.  For  example,  if  the  initial  condition  is  based  on  an  initial  heading  error  in 
the  missile,  then  both  models  yield  the  same  miss  distance  results.  However,  if  the  initial 
condition  stems  from  a  step  in  target  displacement  (this  condition  is  used  to  study  missile 
performance  against  multiple  targets),  then  the  two  block  diagram  topologies,  as  currently 
employed  in  the  open  literature,  will  lead  to  miss  distance  results  that  differ  significantly. 
Consequently,  applying  the  adjoint  method  to  these  two  systems  will  also  give  different 
results.  To  correct  this  inconsistency,  the  necessary  mathematics  is  derived  to  bring  the  two 
formulations  into  alignment  for  any  of  the  initial  state  conditions.  This  then  leads  to 
consistency  in  the  associated  adjoint  models  regardless  of  which  linear  formulation  is  pursued 
by  the  analyst.  In  addition,  closed  form  solutions  of  the  miss  distance  variation  in  terms  of 
flight  time  have  been  derived  for  the  case  when  the  missile  guidance  loop  may  be 
approximated  dynamically  by  a  first  order  lag  term.  Finally,  simulation  results  and  design 
curves  are  generated  for  more  realistic  guidance  systems. 
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2.  Non-linear  Engagement  Equations 

2.1  Planar  Engagement  Model 

Consider  a  missile-target  engagement  under  the  following  assumptions;  the  engagement  is 
confined  to  the  X-Y  plane,  the  force  due  to  gravity  is  ignored,  and  the  missile  velocity  and  the 
target  velocity  remain  constant.  The  geometry  of  the  engagement  is  depicted  in  Figure  1. 

The  non-linear  differential  equations  describing  the  motion  of  the  target  are  [1], 


XT  =  -VT  cos (J3) 

(1) 

Y  =  VT  sin (/?) 

(2) 

^  n 

P=  — 

V 

T 

(3) 

where  (XT ,  Yr  )  defines  the  target  position,  VT  is  the  constant  target  speed,  nT  is  the  target 
acceleration  while  ft  refers  to  the  target  flight  path  angle  as  shown  in  Figure  1.  The  dots  over 
the  variables  denote  differentiation  with  time.  The  initial  conditions  are  given 
by  (0)  =  XT0 ,  Yt  (0)  =  Yto  and  J3( 0)  =  ft, . 


Figure  1.  Engagement  Geometry 
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Similarly,  the  non-linear  differential  equations  describing  the  motion  of  the  missile  are. 


X  =v 

M  MX 

(4) 

Y  =V 

M  MY 

(5) 

=  -n  sin(/l) 

(6) 

=  n  cos(/l) 

(7) 

Here,  (xm,ym)  defines  the  instantaneous  position  of  the  missile  while  ( VMX , VMY ) are  the 
components  of  its  velocity  vector.  Furthermore,  nc  denotes  the  missile  commanded 
acceleration  and  A  is  the  line  of  sight  angle  as  shown  in  Figure  1.  The  initial  conditions  are 
given  by  xM(0)  =  xMO,  fm(0)  =  fmo,  vmx(0)  =  vucos (l  +  he  +  a.)  and  vMY(0)  =  vusm (L  +  he  +  A)  where 
VM  denotes  the  constant  missile  speed,  L  is  the  missile  lead  angle  associated  with  the  collision 
triangle  and  HE  refers  to  the  initial  deviation  of  the  missile  from  the  collision  triangle 
commonly  known  as  the  heading  error.  The  lead  angle  L  can  be  found  by  application  of  the 
law  of  sines  on  the  collision  triangle  yielding  the  formula 

L  =  sin {VT  sin(/?  +  A)  IV M} .  (8) 

To  find  the  missile  acceleration  components,  it  is  necessary  to  determine  the  components  of 
the  relative  missile-target  separation.  Let  the  components  of  the  relative  missile-target 
separation  be  expressed  as 


K=xt-xm, 

A  =Y  -  Y  . 

Y  T  M 

Consequently,  the  range  R  between  missile  and  target  is  given  by 

R  '  (  V  A  )  . 


(9) 

(10) 


(11) 


while  the  line  of  sight  angle  A  is  given  by 


A  =  tan  1 


f 


(12) 


The  line  of  sight  rate  can  be  easily  derived  from  this  expression  by  differentiating  with  respect 
to  time,  giving 
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A  = 


A  A  -A  A 

X  Y  X 


Y 


R2 


(13) 


The  closing  velocity  is  defined  as  the  negative  rate  of  change  of  the  missile  target  separation, 
that  is,  Vc  =  -R ,  which,  after  using  eq  (11)  above,  leads  to  the  relation. 


y  =  -(AA+AA)  (14) 

'  R 

Consequently,  the  magnitude  of  the  missile  guidance  command  jq  can  be  found  from  the 
definition  for  proportional  navigation  guidance, 

n  =  N'V X  ,  (15) 

C  c  y 


where  N  is  a  given  constant. 

If  we  model  the  actual  acceleration  of  the  missile  n  by  a  first  order  lag  term,  then 


n  1 

_ z  _  _ 

nc  1  +  is 


(16) 


A  2D  simulation  model  is  easily  developed  using  the  above  equations.  For  this  study, 
MATLAB / Simulink  was  used  to  develop  the  simulation  model.  The  top  level  Simulink  model 
is  presented  in  Figure  2.  The  details  of  each  subsystem  of  the  model  are  given  in  Appendix  A. 
Figure  3  presents  a  screen  shot  of  the  Graphical  User  Interface  (GUI)  used  in  this  study.  The 
GUI  was  specifically  designed  as  an  aid  in  running  the  simulation.  The  MATLAB  software 
underpinning  the  construction  of  the  GUI  is  also  presented  in  the  Appendix. 

At  the  GUI  level,  the  missile  and  target  parameters  may  be  entered  by  the  analyst.  On  the  left 
of  the  GUI,  the  analyst  can  enter  initial  target  position,  target  speed  and  flight  path  angle.  On 
the  right  of  the  GUI,  the  analyst  can  enter  initial  missile  position,  missile  speed,  effective 
navigation  constant  and  guidance  loop  time  constant.  Inputs  for  the  simulation  include  target 
manoeuvre  (NT),  error  in  the  initial  heading  angle  (FIE)  and  jump  in  the  target  position 
(Displace)  at  some  time  (THOM)  prior  to  intercept. 
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Figure  2  Two  Dimensional  Engagement  Model  in  Simulink 


Figure  3.  Graphical  User  Interface  to  run  the  2D  Engagement  Model 
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2.2  Simulation  Results 

The  Simulink  model  given  above  may  now  be  used  to  conduct  simulation  studies  of  the 
missile/  target  engagement  problem.  The  simulation  studies  will  be  based  on  the  following 
three  simulation  inputs; 

(a)  a  step  in  target  manoeuvre 

(b)  an  error  in  initial  heading  angle 

(c)  a  step  in  target  position 

Interest  in  (c)  above  stems  from  the  analysis  of  multi-target  scenarios  as  is  described  in 
Section  2.2.3. 


2.2.1  Step  in  target  manoeuvre 

Consider  the  case  in  which  the  only  disturbance  is  a  3 g  target  manoeuvre  starting  at  time 
f  =  0 .  The  nominal  values  for  target  and  missile  parameters  may  be  easily  identified  in  the 

GUI  given  in  Figure  3.  In  this  scenario,  the  missile  and  target  are  initially  on  a  collision  course 
and  flying  along  the  downrange  component  of  the  earth  fixed  co-ordinate  system.  Thus,  the 
target  velocity  vector  is  initially  along  the  line  of  sight  and,  at  first,  all  3 g  of  the  target 
acceleration  are  perpendicular  to  the  line  of  sight.  However,  as  the  target  manoeuvres,  the 
magnitude  of  the  target  acceleration  perpendicular  to  the  line  of  sight  reduces  due  to  the 
turning  of  the  target. 

Sample  missile/  target  trajectories  for  this  case  with  effective  navigation  ratios  of  3  and  5  are 
shown  in  Figure  4.  It  is  clear  from  the  figure  that  the  higher  effective  navigation  ratio  causes 
the  missile  to  lead  the  target  slightly  more  than  the  lower  navigation  ratio  case. 

Figure  5  presents  the  respective  missile  acceleration  profiles  obtained  from  the  simulations. 
Note  that,  although  both  acceleration  profiles  are  monotonically  increasing  for  most  of  the 
flight,  the  higher  effective  navigation  ratio  case  leads  to  less  acceleration  requirement  of  the 
missile  towards  the  end  of  flight.  Also  noteworthy  from  the  plot  is  the  observation  that  the 
peak  acceleration  required  by  the  missile  to  hit  the  target  is  significantly  higher  than  the 
manoeuvre  level  of  the  target  (3g). 

2.2.2  Heading  error 

Next,  consider  the  case  in  which  the  only  disturbance  is  a  20  degree  error  in  the  initial  heading 
angle,  that  is,  HE  =  -20  deg.  Again  the  simulation  was  run  for  two  values  of  the  effective 
navigation  ratio.  Sample  trajectories  for  effective  navigation  ratios  of  3  and  5  are  presented  in 
Figure  6.  From  the  figure,  it  is  apparent  that  initially  the  missile  is  flying  in  the  wrong 
direction  because  of  the  heading  error.  Gradually  the  guidance  law  forces  the  missile  to  head 
towards  the  target.  It  is  interesting  to  note  from  the  figure  that  the  larger  effective  navigation 
ratio  enables  the  missile  to  remove  the  initial  heading  error  more  rapidly,  thus  leading  to  a 
tighter  missile  trajectory  in  this  case. 
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In  Figure  7  is  plotted  the  resultant  missile  acceleration  profiles  for  each  case.  From  the  figure, 
it  is  observed  that  the  faster  removal  of  heading  error  in  the  higher  effective  navigation  ratio 
case  is  associated  with  larger  missile  accelerations  at  the  beginning  of  flight 
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Figure  4:  Missile  performance  against  manoeuvring  target 
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Figure  6:  Missile  performance  when  initial  heading  angle  is  in  error 
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Figure  7:  Missile  performance  when  initial  heading  angle  is  in  error 
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2.2.3  Step  in  target  position 

The  non-linear  engagement  simulation  program  may  be  used  to  analyse  the  missile  response 
to  a  sudden  jump  in  target  position  at  a  certain  time  prior  to  intercept.  This  is  useful  in 
supporting  simulation  studies  concerning  the  effects  of  seeker  resolution  on  the  performance 
of  the  missile  when  confronted  with  multiple  targets. 

Consider  the  scenario  shown  in  Figure  8  [1],  On  the  left  is  the  missile  engaging  two  targets 
which  are  flying  in  formation.  Both  targets  have  the  same  speed  and  are  separated  by  a 
spacing  of  Ad  metres.  It  is  assumed  here  that  the  power  centroid  is  located  half  way  between 
the  targets.  In  this  case,  the  input  parameter,  "Displace"  in  the  GUI  of  Figure  3  is  equivalent 
to  Ad  /  2  metres.  For  the  simulation,  it  is  assumed  that  the  missile  is  initially  on  a  collision 
course  with  the  power  centroid.  At  a  certain  time  to  go  before  intercept  with  the  power 
centroid,  seeker  resolution  occurs  and  the  missile  is  presented  with  the  true  target,  that  is. 
Target  1.  At  this  time,  it  will  appear  to  the  missile  as  if  the  target  position  jumps  from  the 
power  centroid  to  Target  1.  From  a  simulation  perspective,  it  is  only  necessary  to  model  the 
target  currently  seen  by  the  missile.  Therefore,  for  most  of  the  flight,  the  missile  will  be 
guiding  on  the  power  centroid  and  for  the  rest  of  the  flight,  following  seeker  resolution,  the 
missile  will  be  guiding  on  Target  1.  This  is  represented  in  the  simulation  by  the  target  position 
being  updated  by  a  step  in  target  displacement  at  a  given  time  to  go  prior  to  intercept 
(THOM). 

The  non-linear  2D  engagement  model  is  used  to  generate  sample  simulation  results  of  this 
type  of  engagement.  Relevant  simulation  results  are  presented  in  Figure  9  and  10  below. 
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Figure  8:  Missile  engaging  two  targets  flying  information 
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Figure  9:  Missile  performance  in  multi-target  case 
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3.  Linearised  Engagement  Equations 

In  this  section,  analytical  tools  are  employed  to  gain  further  insight  into  the  performance  of 
the  interceptor  under  different  engagement  conditions  and  to  better  understand  its  dynamic 
response.  In  particular,  the  adjoint  method  is  applied  to  a  simplified  version  of  the 
engagement  equations.  The  simplified  equations  of  motion  are  derived  using  linearization 
about  the  line  of  sight  angle  between  missile  and  target.  The  procedure  follows  closely  that 
outlined  in  Zarchan  [1]. 

First,  define  the  relative  separation  between  the  missile  and  target  perpendicular  to  the  fixed 
reference  as  defined  in  Figure  11, 


y=Y -Y  . 

STM 


The  relative  acceleration  is  expressed,  by  inspection  of  the  figure,  as 

y  =  nr  cos  j3  -  nc  cos  X . 


(17) 


(18) 


Figure  11:  Engagement  Geometry  for  Linearisation 

For  small  flight  path  angles,  that  is,  near  head-on  or  tail  chase  case,  the  cosine  terms  are 
approximately  unity,  and  the  previous  equation  (18),  reduces  to 

y  =  nT-nc.  (19) 
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Similarly,  the  line  of  sight  angle,  which  is  given  by, 

•  ^  y 

sin  A,  =  — , 
R 


(20) 


may  also  be  linearised,  using  the  small  angle  approximation,  yielding  the  expression 


For  a  head-on  case,  the  closing  velocity  is  approximated  by 

V  =  V  +  V 

and  in  a  tail  chase  situation,  the  closing  velocity  is  approximated  by 

V  =V  -V  . 

C  M  T 


(21) 


(22) 


(23) 


Consequently,  for  constant  missile  and  target  speed,  the  closing  velocity  may  be  treated  as  a 
positive  constant.  However,  the  closing  velocity  has  previously  been  defined  as  the  negative 
derivative  of  the  range  from  missile  to  target,  that  is,  Vc  =  -R ,  and  since  the  range  must  go  to 

zero  at  the  end  of  flight,  we  can  approximate  the  range  equation  with  the  time  varying 
relationship. 


R  =  Vc(tp-t).  (24) 

In  the  above  expression,  t  denotes  the  current  time  and  tF  is  the  total  flight  time  of  the 
engagement.  Note  that  tF  is  a  constant  here. 

The  linearised  miss  distance  is  defined  to  be  the  relative  separation  between  missile  and  target 
at  the  end  of  flight,  namely. 


MD  =  y(tF). 


(25) 


Using  equations  (19),  (21),  (24)  and  (25),  and  adding  the  information  given  by  equations  (15) 
and  (16),  we  are  able  to  build  a  block  diagram  model  for  the  linearised  version  of  the  homing 
equations.  This  model  is  displayed  in  Figure  12  below.  In  the  diagram,  the  symbol  s  refers  to 
differentiation  in  the  frequency  domain  using  Laplace  Transform  terminology. 

In  order  to  implement  this  block  diagram  model  in  Simulink,  it  is  prudent  to  combine  the 
single  s  block  with  the  transfer  function  block  given  by  1  /(I  +  is)  •  This  is  possible  in  this  case 
as  N'V  is  a  constant.  Thus,  following  this  process,  the  revised  block  diagram  model  for  the 
linearised  equations  is  presented  in  Figure  13.  This  model  is  designated  here  as  Model  A.  Note 
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that  the  block  diagram  displayed  in  Figure  13  now  contains  an  inner  closed  loop  block 
structure.  This  alternative  block  structure  is  mathematically  equivalent  to  the  single  block 
with  transfer  function  s/(l  +  zy).  The  process  is  known  as  reducing  the  block  to  its 
fundamental  closed  loop  form  and  has  been  achieved  here  using  standard  block  diagram 
algebra.  This  alternative  representation  is  important  in  the  sequel. 


MD  =  y(t„) 


Figure  12:  Linearised  Geometry  Model 


MD  =  y{tp) 


Figure  13:  Linear  Model  A 


Model  A  is  typically  employed  in  support  of  analytical  studies  of  the  homing  loop  problem  in 
the  literature  [1,2,3].  This  model  is  also  used  as  the  basis  for  adjoint  analysis  of  the  guidance 
loop.  It  is  shown  in  the  sequel  that  while  this  model  is  useful  for  performance  studies 
associated  with  target  manoeuvre  and  heading  error  effects,  caution  needs  to  be  exercised 
when  using  the  model  as  a  basis  for  analysis  of  multi-target  problems. 

An  alternative  block  diagram  structure  typically  employed  in  the  literature  for  homing  loop 
studies  [1,4,5,6]  may  be  derived  using  the  analytical  expression  for  the  rate  of  change  of  line  of 
sight  angle.  For  the  non-linear  model,  the  rate  of  change  of  the  line  of  sight  angle  is  given  in 
equation  (13).  After  applying  the  linearisation  conditions  stipulated  above,  the  expression  for 
the  line  of  sight  rate  becomes 
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x = _ z _ + _ z _  (26) 

K(tF-t)  vc(tF  -ty  ’ 

or  equivalently, 

^  =  itj  +  y)  (27) 

T  T  .  7  ? 


where  the  term  /  denotes  the  time  to  go  and  is  defined  as  /  =  /  —  / . 

go  go  F 

After  incorporating  the  analytical  expression  for  the  line  of  sight  rate  into  the  linearised  set  of 
equations,  an  alternative  block  diagram  for  the  homing  loop  dynamics  may  be  derived.  This 
alternative  block  diagram,  designated  as  Model  B,  is  displayed  in  Figure  14. 


Figure  14:  Linear  Model  B 


In  the  literature,  these  two  models  (Model  A  and  Model  B)  are  often  used  interchangeably  as 
the  basis  for  a  preliminary  analysis  of  the  performance  of  the  generic  guided  missile  homing 
loop  [1-6].  The  next  section  covers  a  comparison  of  the  two  models  through  the  use  of 
simulation  and  highlights  input  conditions  under  which  the  models  yield  different  results.  For 
these  input  conditions,  caution  should  be  exercised  by  the  analyst  when  adopting  these 
models,  particularly  when  considering  Model  A. 

4.  Comparison  of  the  Linearised  Models 

The  two  linear  models  have  been  implemented  in  Simulink  in  preparation  for  a  comparative 
simulation  study.  The  Simulink  implementation  of  Model  A  and  B  are  presented  in  Figures  15 
and  16,  respectively.  For  the  present  work,  the  simulations  of  the  linear  models  are  carried  out 
under  the  same  input  conditions  considered  in  Section  2.2  above.  This  allows  a  comparison  of 
the  simulation  results  generated  by  the  linear  models  against  those  obtained  by  the  non-linear 
model.  To  do  this,  the  time  of  intercept  displayed  in  the  GUI  after  running  the  non-linear 
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model  is  recorded  as  the  approximate  flight  time,  tF ,  for  the  engagement.  This  flight  time  is 
then  used  in  the  linear  models.  Additionally,  the  same  numerical  values  for  the  parameters 
are  used,  as  shown  in  the  GUI  in  Figure  3. 


Integ2 


Figure  15:  Simulink  Implementation  of  Linear  Model  A 


Figure  16:  Simulink  Implementation  of  Linear  Model  B 

4.1.1  Step  in  target  manoeuvre 

A  step  in  target  acceleration  of  magnitude  3 g  is  applied  to  both  linear  models  and  all  other 
inputs  are  set  to  zero.  The  simulation  results  are  presented  in  Figures  17  and  18  below. 
Figure  17  shows  the  time  history  of  the  relative  displacement  profile  while  Figure  18  shows 
the  achieved  missile  acceleration  profile.  It  is  clear  from  the  figures  that  both  Model  A  and 
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Model  B  simulation  results  agree.  Furthermore,  the  linear  models  provide  a  good 
approximation  to  the  non-linear  engagement  model  in  this  case. 


Figure  17:  Comparison  of  relative  displacement  profile  in  the  case  of  a  3g  step  in  target  manoeuvre 


Time  (s) 


Figure  18:  Comparison  of  missile  acceleration  profile  in  the  case  of  a  3g  step  in  target  manoeuvre 
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4.1.2  Heading  error 

In  this  case,  an  error  in  initial  heading  angle  of  -20  degrees  is  considered  with  all  other  inputs 
set  to  zero.  Again,  the  simulation  results  from  the  linear  models  are  compared  with  those 
generated  by  the  non-linear  engagement  model.  The  heading  error  is  introduced  in  the  linear 
models  as  an  appropriate  initial  condition  on  the  first  integrator  in  the  Simulink  model  of 
Figure  15  and  16,  respectively.  The  results  of  the  simulations  are  presented  in  Figure  19  and  20 
below. 

The  figures  indicate  that  both  linear  models  agree  in  this  case  and  that  they  represent  a  good 
approximation  to  the  dynamics  of  the  non-linear  engagement  model. 


Figure  19:  Comparison  of  relative  displacement  profile  in  the  case  of  a  heading  error  of -20  deg 

4.1.3  Step  in  target  position 

Here,  interest  lies  in  the  response  of  the  missile  to  a  sudden  step  in  target  position  at  a  certain 
time  to  go  prior  to  intercept.  This  scenario  simulates  the  effect  of  target  resolution  by  the  on¬ 
board  seeker  in  a  multi-target  scenario.  For  the  linearised  models,  this  condition  translates 
into  an  appropriate  initial  condition  placed  on  the  second  integrator  (Integl)  in  Figures  15  and 
16. 

Simulink  simulations  were  carried  out  for  both  linear  models  and  then  compared  to  the 
corresponding  non-linear  results,  as  shown  in  Figures  20  and  21  below.  These  plots  show  that 
Model  A  results  do  not  agree  with  Model  B  results.  Furthermore,  the  plots  demonstrate  that 
Model  B  results  are  a  good  approximation  to  the  simulation  results  generated  by  the  non¬ 
linear  model.  However,  the  same  cannot  be  said  of  Model  A  in  this  case. 
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Figure  20:  Comparison  of  missile  acceleration  profile  in  the  case  of  a  heading  error  of -20  deg 

Consequently,  the  observations  gleaned  from  the  foregoing  simulations  suggest  that  Model  B 
is  a  sufficiently  accurate  linear  model  when  considering  the  missile  target  engagement 
problem.  And  this  is  so  regardless  of  the  three  input  conditions  considered.  However,  caution 
needs  to  be  exercised  when  using  Model  A  for  this  purpose  as  it  has  been  shown  via 
simulation  to  produce  erroneous  results  under  one  of  the  stipulated  input  conditions  (see 
Figure  21  and  22). 
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Figure  21:  Comparison  of  relative  displacement  profile  in  the  case  of  a  jump  in  target  position  of 
60  m 
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Figure  22:  Comparison  of  missile  acceleration  profile  in  the  case  of  a  jump  in  target  position  of 
60  m 
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5.  Derivation  of  a  Correction  for  Model  A 


In  this  section,  an  analysis  is  carried  out  to  explore  whether  Model  A  can  be  brought  into 
alignment  with  Model  B  when  the  input  condition  is  a  step  in  target  displacement.  Consider 
the  block  diagram  of  Model  A  given  in  Figure  13  and  reproduced  in  Figure  23  for 
convenience.  In  Figure  23,  the  states  of  the  system  have  been  included  in  the  block  diagram, 
designated  by  y ,  as  has  the  desired  input  condition  yic . 


MD  =  y(tF) 


Figure  23:  Linear  Model  A  with  system  states  displayed 

By  inspection  of  the  diagram,  the  state  space  equations  may  be  easily  written  down.  They  are. 


i,  =  nT  -  N’Vcx}  , 

(28) 

x2  =  x,  , 

(29) 

K  x,  , 

*3=  ( Tr  ;  X3). 

(30) 

In  this  case,  the  initial  conditions  are,  x{  (0)  =  0/  X2  (0)  =  y.  and  x3  (0)  =  0  •  Thus,  for  time 
zero,  equation  (30)  yields 


*M  = 


i  yk 

T  VctF  ' 


(31) 


Consequently,  the  initial  value  of  the  relative  acceleration  j;(0)  will  have  the  magnitude 
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y(0)  =  xl  (0)  =  nT 


N'vic 


(32) 


This  is  physically  impossible  as  it  implies  a  jump  in  missile  achieved  acceleration  whenever 
the  initial  condition  on  y  is  non-zero.  To  counteract  this  effect,  a  requirement  is  placed  on  the 
filter  initial  condition  as  follows. 


x3(0)  = 


y* 

vctF ' 


(33) 


When  this  correction,  in  the  form  of  a  non-zero  initial  condition  on  the  filter,  is  applied  to 
Model  A,  the  simulation  results  generated  by  both  linear  models  agree  as  shown  in  Figure  24 
and  25  below. 

Thus,  in  conclusion,  any  of  the  linear  models.  Model  A  or  Model  B,  derived  previously  may  be 
used  for  linear  analysis  of  the  missile  target  engagement  problem  provided  that  Model  A  is 
used  in  conjunction  with  the  correction  factor  derived  above.  It  is  noteworthy  to  point  out  that 
in  the  case  when  yjc  is  zero,  no  initial  correction  of  the  filter  is  necessary. 


Figure  24:  The  relative  displacement  profile  generated  by  the  linear  models  now  agree  for  the 
initial  displacement  case 
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Time  (s) 


Figure  25:  The  missile  acceleration  profile  generated  by  the  linear  models  now  agree  for  the  initial 
displacement  case 


6.  Analytical  Formulas 

The  linearisation  of  the  engagement  model  is  important  for  two  reasons.  Firstly,  with  a  linear 
model,  powerful  computerised  techniques,  such  as  the  adjoint  method,  may  be  utilised  to 
analyse  the  missile  guidance  system  both  statistically  and  deterministically  in  one  computer 
simulation.  Moreover,  with  this  technique,  error  budgets  are  automatically  generated  so  that 
key  system  drivers  may  be  identified  and  a  balanced  system  design  can  be  achieved. 
Secondly,  under  special  circumstances,  the  linear  engagement  model  is  mathematically 
amenable  to  analytical  solutions.  These  solutions  can  assist,  during  the  preliminary  phases  of 
a  missile  design,  in  gaining  insights  for  system  sizing.  Furthermore,  the  form  of  the  analytical 
solutions  will  provide  clues  on  how  key  parameters  may  influence  system  performance. 

In  this  section,  closed-form  solutions  are  derived  for  the  three  important  cases  that  were 
considered  above  for  the  engagement  simulations.  The  aim  is  to  derive  closed-form  solutions 
associated  with  Model  A  and  Model  B  independently  and  then  to  show  that  these  solutions 
are  mathematically  equivalent  provided  that  the  filter  initial  condition  is  included  in  Model  A. 
It  has  been  shown  in  the  literature  [1]  that  closed  form  solutions  for  the  linear  engagement 
model  may  be  obtained  more  easily  using  the  adjoint  method.  Therefore,  this  is  the  method 
employed  in  the  derivation  process. 
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6.1  Closed-Form  Solutions  for  Model  A 

The  adjoint  block  diagram  of  Model  A  is  presented  in  Figure  26.  This  model  was  constructed 
from  the  forward  loop  model  given  in  Figure  13  using  the  usual  rules  of  adjoint  construction 
given  in  references  [1,  3].  Note  that  x  in  the  diagram  refers  to  the  filter  initial  condition 

derived  in  the  last  section  and  is  expressed,  in  this  case,  as 


*  =  (34) 

*  vctF' 

Also  in  the  diagram,  the  initial  condition  for  the  heading  angle  error  is  defined  as 


HA  =  -V  HE. 

M 


(35) 


MD 


NT 


Figure  26: 


Adjoint  Model  A 


In  this  form,  the  state  equations  associated  with  the  adjoint  block  diagram  (Figure  26)  are 
mathematically  challenging  for  analytical  work.  However,  using  block  diagram  algebra  [1], 
the  structure  of  this  block  diagram  may  be  conveniently  reduced  to  a  form  such  as  to  allow 
easy  mathematical  treatment  of  the  problem.  The  reduced  form  of  the  block  diagram  is 
presented  in  Figure  27  below.  In  the  diagram,  h(t)  is  designated  as  an  adjoint  output  signal 
while  the  weighting  function  w(t),  after  frequency  domain  transformation,  is  defined  as 

N'  (36) 

5(1  +  ZS ) 
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Figure  27:  Reduced  Form  of  Adjoint  Model  A 


Using  the  convolution  integral,  the  adjoint  output  is  related  to  the  input  by 

t 

h(t )  =  j*  w(ji)[8{t  -rj)-h(t  ~rj)]drj . 

0 

After  taking  the  Laplace  Transform  of  the  preceding  equation,  one  obtains 

-dA  =  W(s)[\-H{s)}. 

as 

However,  since 

_dH=J[l  H(s)]'  (39) 

els  as 

the  technique  of  separation  of  variables  is  invoked  to  arrive  at 

J  ALAI  =  ]>(.)*  .  <«> 

The  solution  to  equation  (40)  is  given  by 

l-tf(s)  =  Cexp(JlL(s>/s)  ,  (41) 

where  C  is  a  constant  of  integration.  This  constant  can  be  evaluated  by  first  noting  from 
Figure  27  that  the  miss  distance  due  to  a  step  in  target  displacement  of  magnitude  y  is  given 

by  (using  frequency  domain  notation) 


(37) 


(38) 
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L{MDv}  =  ^[\-H{s)\  .  (42) 

y*  s 

In  the  above  expression,  the  term  l{  }  refers  to  the  Laplace  transform  of  its  argument.  Let  the 
magnitude  of  the  step  be  unity.  Then,  in  the  time  domain,  the  miss  distance  due  to  a  unit  step 
in  target  displacement  at  time  zero  will  be  unity.  Thus,  invoking  the  initial  value  theorem 
leads  to 


MDJ  0)  =  1  = 


Inn 

5  -»  00 


V  ^ 


Therefore,  C  is  chosen  to  satisfy  the  following  relation, 

51™(X)Cexp(JlL(5V5)  =  l  . 


(43) 


(44) 


In  this  case,  since  JV(s)  is  given  by  equation  (36),  then  it  follows  that  C  =  1,  and  equation  (41) 
is  reduced  to 


H(s )  =  exp(J  W(s)ds)  . 


(45) 


Consequently, 


fw(s)ds  =  N’  J 


ds 

5(1  +  IS) 


N 


i 


5  (5  +  1/r) 


]  ds  , 


f 


=  N'  In 


A 


s  +  l/r 


(46) 

(47) 

(48) 


Now,  after  substitution  of  this  expression  into  equation  (45)  above,  one  obtains  the 
fundamental  relationship 


1  ~H(s) 


VA  +  l/r  ) 


(49) 


Consequently,  from  inspection  of  the  adjoint  block  diagram  given  in  Figure  27  above,  it  is  easy 
to  write  down  expressions  (in  frequency  domain  notation)  for  each  of  the  desired  outputs  as 
follows.  The  miss  distance  due  to  a  step  manoeuvre  of  magnitude  fj  is  given  by 
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l[MDnt]_\-H{s)_  1 


^  s  T 


nT  s 3  s^s  +  l/r. 


The  miss  distance  due  to  an  initial  heading  error  is  given  by 


VMD^  v  [1  -H(S)\_-VM 
1  HE  1  M  s 2  s2 


^  s  ^ 


s  +  l/r 


(50) 


(51) 


For  a  step  in  target  displacement,  the  miss  distance  formula  is  composed  of  two  parts,  namely, 
that  due  to  y  and  that  due  to  x  ■  The  part  due  to  y  is 

y  ic  ic  -L  y  ic 


l{MD^}  =  =  U  s  " 

yic  s  s{s  +  1/t) 


(52) 


while  the  part  due  to  x  is,  according  to  Figure  27, 


l<MD  _[\-H(sW(s)  _N 

VrX.  S  TS 

C  ic 


'  r 


Vs  +  l/rJ 


1 

S  +  H  T 


(53) 


Therefore,  for  an  effective  navigation  ratio  of  3,  the  inverse  Laplace  Transform  of  the  above 
expressions  is  taken,  yielding  the  analytical  expressions  in  equations  (54)  and  (55).  Note  that 
the  time  variable  appearing  in  the  expressions  is  adjoint  time  which  may  be  interpreted  as 
total  flight  time  £  for  the  engagement  problem. 


MD, 


n „ 


0.5  t]e 


and 


MD 


HE 

HE 


-V..tFe 

M  F 


-tF  /  T 


6  t  \ 
1-^ 

V  2r ) 


(54) 


(55) 


Furthermore,  the  miss  distance  due  to  a  step  in  target  displacement  may  be  determined  by  the 
addition  of  the  inverse  Laplace  Transform  of  equations  (52)  and  (53).  The  closed  form 
expression  for  this  is 


MD 

_ 

l  y,c  ) 


Tot  |V  =3 


MD 

, _ y 

y. 


+- 


MD 

xi 

yic 


(56) 
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where 

MD 

yic 

y  &  N=  3 

and 


(58) 


In  deriving  equation  (58)  above,  use  was  made  of  equation  (34).  Consequently,  adding  the 
above  expressions  yields  the  miss  distance  due  to  a  step  in  target  displacement  of  magnitude 
y  ,  namely 


=  [1.5—  -0.5  *y]e~ 

T  T~ 


2  — +  0.5— ]A 

T  T 


[1-0.5— ]e“',/r  . 

T 


(59) 


6.1.1  Adjoint  of  Model  A  in  Simulink 

A  Simulink  model  of  the  adjoint  system  of  Model  A  (Figure  26)  is  shown  in  Figure  28.  This 
model  is  used  to  generate  simulation  data  for  comparison  with  the  closed  form  solutions 
given  previously.  The  comparisons  are  displayed  in  the  following  figures. 


Figure  28:  Adjoint  Model  A  in  Simulink 
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Figure  29:  Adjoint  Simulation  and  Formula  agree  for  MDnt 


Figure  30:  Adjoint  Simulation  and  Formula  agree  for  MDh  e 
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Figure  31:  Adjoint  Simulation  and  Formula  agree  for  MDyic 

The  plots  confirm  that  the  closed  form  solutions  and  the  Simulink  simulation  results  based  on 
the  linear  adjoint  of  Model  A  agree  for  the  particular  cases  considered  here.  Thus,  either 
approach  may  be  used  for  preliminary  system  design.  Recall  that  the  closed  form  solutions  are 
based  on  a  single  time  constant  guidance  loop.  However,  for  higher  order  systems,  the 
equations  become  mathematically  intractable.  Thus,  using  the  simulation  model  allows  the 
analyst  added  flexibility  for  studying  more  complex  homing  loop  systems  with  minimum 
extra  effort  due  to  the  block  diagram  features  of  Simulink.  This  aspect  shall  be  considered 
further  in  a  later  section  of  the  report. 

6.1.2  Comparison  with  non-linear  results 

Here,  the  non-linear  engagement  model,  as  presented  in  Section  1,  is  used  in  a  multi-run  mode 
to  generate  the  miss  distance  profile  as  a  function  of  flight  time  or  time  to  go.  This  will  enable 
a  comparison  of  the  non-linear  simulation  results  with  those  obtained  using  either  the  linear 
adjoint  model  in  Simulink  or  the  closed  form  solutions  derived  earlier.  For  the  non-linear 
model,  the  approximate  miss  distance  calculation  method,  as  employed  by  Zarchan  [1],  is 
utilised.  Furthermore,  in  order  to  obtain  the  miss  distance  profile  as  a  function  of  flight  time 
using  the  non-linear  simulation,  a  range  to  go  parameter  is  defined  and  utilised  in  the  process. 

The  miss  distance  generated  using  the  non-linear  model  for  a  step  in  target  manoeuvre  of 
magnitude  3g  is  shown  in  Figure  32  for  different  flight  times.  Superimposed  on  this  plot  is  the 
miss  distance  profile  computed  from  the  closed  form  solution  for  this  scenario.  It  is  clear  from 
this  plot  that  the  closed  form  solution  provides  sufficiently  accurate  results  for  low  values  of 
flight  time.  As  the  flight  time  increases,  the  plot  shows  that  the  non-linear  results  tend  to 
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diverge.  According  to  linear  theory,  the  miss  distance  tends  to  zero  for  flight  times  that  are 
approximately  greater  than  ten  times  the  missile  response  time  constant.  This  is  clearly 
evident  in  the  figure  as  the  missile  time  constant  is  0.3  in  this  case.  Consequently,  beyond  a 
flight  time  of  3  seconds,  linear  theory  is  less  accurate  in  this  case. 

In  Figure  33,  the  non-linear  simulation  results  are  compared  with  those  obtained  using  the 
closed  form  solution  for  the  case  of  an  error  in  initial  missile  heading  angle.  In  this  case,  good 
agreement  is  obtained  for  most  of  the  flight  times  considered.  Large  discrepancies  occur  when 
the  total  flight  time  is  very  small.  This  stems  from  the  fact  that  under  such  conditions,  the 
engagement  scenario  is  significantly  non-linear.  As  a  final  comparison  with  the  non-linear 
simulation,  the  case  of  a  step  in  target  displacement  is  considered.  In  this  case,  the  non-linear 
simulation  results  were  obtained  by  looping  over  different  values  of  homing  time,  or  THOM 
in  Figure  3.  Again,  the  closed  form  solutions,  based  on  the  linear  analysis,  agree  reasonably 
well  with  the  non-linear  simulation  results  generated  by  the  Simulink  model  (see  Figure  34). 
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Figure  32:  Formula  agrees  closely  with  non-linear  simulation  for  MDNT 
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6.2  Closed-Form  Solutions  for  Model  B 

In  this  section,  linear  Model  B  is  investigated  using  analytical  means.  The  forward  loop  block 
diagram  of  Model  B  is  given  in  Figure  14.  Applying  the  adjoint  construction  rules  to  this  block 
diagram  yields  the  adjoint  block  diagram  of  Model  B  as  presented  in  Figure  35. 

In  Figure  35,  the  outputs  of  two  of  the  integrators  in  the  block  diagram  are  designated  as 
Z  and  2,  -  These  are  clearly  indicated  in  the  figure.  From  inspection  of  the  block  diagram,  the 

following  expressions  may  be  deduced,  namely 

i,  =  z2  +  tz2 ,  z2  (0)  =  1 .  (60) 

Consequently,  after  integrating  this  equation,  one  obtains  the  expression 

21  =  tZ2  +  >  (61) 

where  C,  is  a  constant.  The  value  of  this  constant  can  be  determined  by  noting  from  Figure  35 

that  z  (t)  is  directly  linked  to  the  miss  distance  due  to  a  heading  angle  error.  Recall  that  the 

time  variable  in  this  case  (adjoint  model)  is  interpreted  as  flight  time.  Thus,  for  zero  flight 
time,  the  following  must  be  true,  namely,  z  (0)  =  C,  =  0  •  Consequently,  equation  (61)  reduces 

to 

z,  =  tz2  .  (62) 


Figure  35:  Adjoint  Model  B 

Hence,  making  use  of  equation  (62)  above,  and  after  introducing  the  expression  for  the 
weighting  function  defined  earlier  in  equation  (25),  the  block  diagram  displayed  in  Figure  35 
may  be  judiciously  manipulated  into  the  block  structure  presented  in  Figure  36. 
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For  ease  of  analysis,  the  block  diagram  displayed  in  Figure  36  is  re-cast  into  the  standard 
input-output  form  shown  in  Figure  37  below.  Note  that  care  must  be  exercised  when  doing 
this  since  the  block  diagram  contains  variables  from  different  domains,  that  is,  from  the  time 
domain  (/)  and  the  frequency  domain  (.s1) .  In  the  block  diagram  of  Figure  37,  two  extra 

variables  have  been  included  as  adjoint  signals  of  interest,  namely,  z  (t)  and  z4  (t)  ■ 


Figure  36:  Equivalent  block  diagram  of  adjoint  Model  B 


Figure  37:  Equivalent  block  diagram  in  standard  form 

From  Figure  37,  the  following  time  domain  expressions  may  be  deduced, 

=  ,  (63) 

t  Jo 
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(64) 


where  f  (t)  is  defined  as 


/(O  -  L~l{sW(s)} 


(65) 


Taking  the  Laplace  Transform  of  equation  (63)  and  (64)  yields  the  following  relationships. 


d2Z4(s) 


ds2 


ZM  = 


=  sW(s)Zl(s)  , 


1  ~ZA{s)  Z3(s) 


(66) 


(67) 


Furthermore,  taking  the  Laplace  Transform  of  equation  (62),  yields  the  relation 

dZ2(s)  (68) 


ZX*)  = 


ds 


After  some  algebra,  the  following  differential  equation  for  Z  (s)  is  derived. 


dZt 

ds 


-  + 


--W(s) 

^  J 


Z  =0  . 


(69) 


This  differential  equation  may  be  readily  solved  using  separation  of  variables  to  yield 

Z,(S)  =  CeN(^  ,  (70) 

S' 


where  C1  is  a  constant.  The  value  of  the  constant  may  be  determined  as  follows.  From  the 

block  diagram  in  Figure  37,  the  expression  for  the  miss  distance  due  to  target  displacement 
may  be  derived.  This  has  the  form 

MDyt  =YT*  z2  (t)  .  (71) 

Therefore,  for  a  unit  step  in  target  displacement  at  zero  flight  time,  the  following  condition 
must  be  true,  namely. 


z2(0)  =  l 


(72) 
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Note  the  consistency  with  equation  (60)  above,  where  the  same  condition  arose  as  a  result  of 
the  impulsive  input. 


After  substituting  the  expression  for  W  (.S’),  as  defined  in  equation  (25),  into  equation  (70) 
above,  and  after  some  algebra,  the  following  expression  is  obtained. 


zAs)  = 


c 

s 2 


s  +  1/t 


(73) 


From  equation  (62) 


Z2  (0  ~ 


Z,(0 


(74) 


Thus  at  zero  flight  time,  the  following  expression  is  valid, 

lim  z,(0 

t  —>  o  t  ' 


(75) 


Using  L'FlopitaTs  rule  and  making  use  of  equation  (72)  yields  the  expression, 

lim  dzt 


1 


dt 


(76) 


Now,  the  Laplace  Transform  of  dzx  /  dt  is  given  by 

l{^4  =  sZ,(s)-2l(  0)  . 
at 


(77) 


However,  as  was  shown  earlier,  z  (0)  =  0  ■  Hence,  the  above  expression  reduces  to 


lyy\  =  sZAs)  ■ 
dt 


(78) 


Consequently,  returning  to  equation  (76)  and  invoking  the  initial  value  theorem  leads  to  the 
following  expression 


,  lim  2  ry  ,  , 

1  =  v  s  Z  (s) 

S  — >  00  1  v  ' 


(79) 


Hence,  after  substituting  equation  (70)  into  equation  (79),  and  making  use  of  equation  (25),  the 
above  expression  reduces  to 
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1=  c, 

S  — »  oo  2 


S  +  H  T 


(80) 


This  implies  that  C2  =  1  •  Consequently,  the  expression  for  Z  (.S')  is 


ZXs)  = 


1 


S' 

s  +  1/t 


(81) 


This  expression  may  now  be  used  to  deduce  the  Laplace  Transform  of  the  miss  distance 
outputs  of  interest  highlighted  in  Figure  37  above,  namely,  the  miss  distance  due  to  target 
manoeuvre. 


n„ 


MDnt  (s1)  =  -J- 


s  +  l/r 


(82) 


and  the  miss  distance  due  to  an  initial  heading  error. 


-V  HF 

MDha  (s)  =  — 


s 

s  +  Mt 


(83) 


Finally,  the  miss  distance  due  to  target  displacement  is  found  by  making  use  of  equation  (74). 
After  time  domain  conversion,  this  may  be  expressed  as 


MDYT(tF)  =  YT*zXtF)!tf 


(84) 


where  t  denotes  flight  time. 


6.2.1  Closed  Form  Solutions  for  N'=3 

Here  closed  form  solutions  are  derived  for  a  particular  value  of  the  navigation  ratio,  namely, 
N  =  3  -  Substitution  of  this  value  into  the  above  expressions,  and  after  taking  the  inverse 
Laplace  Transform,  yields  the  following  closed  form  solutions. 


MD 


n. 


mdha 

HE 


MD 


YT 


YT 


-  =0.5  t;e-hlT  , 

N'= 3 

(85) 

=  -V„t„e-“\  1-L)  , 

2  T 

(86) 

=  e~th  lt  (\  -  • 

(87) 

N'= 3 


36 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2845 


These  agree  exactly  with  the  corresponding  formulas  derived  earlier  and  based  on  Model  A. 
Closed  form  solutions  may  also  be  derived  for  other  values  of  N'  provided  they  are  integers. 
In  fact,  several  of  these  closed  form  solutions  have  been  derived  and  are  reproduced  in 
Appendix  B.  From  the  closed  form  solutions  given  above  (and  in  Appendix  B),  it  is  clear  that 
the  miss  distance  due  to  a  step  in  target  displacement  is  related  to  the  miss  distance  due  to  a 
heading  error.  This  relationship  is  investigated  further  in  Appendix  C. 

However,  for  non-integer  values  of  N',  the  mathematics  is  complex.  In  this  case,  the  best 
approach  for  a  solution  would  be  to  use  simulation.  A  Simulink  representation  of  Model  B  is 
displayed  in  Figure  38  below.  As  the  simulation  results  agree  exactly  with  those  previously 
produced  using  Model  A  (see  Figures  29-31),  only  those  results  generated  using  Model  B  will 
be  presented  here. 


Figure  38:  Adjoint  Model  B  in  Simulink 


The  miss  distance  due  to  a  target  manoeuvre  of  magnitude  3y  (that  is,  output  3  in  Figure  38)  is 
presented  in  Figure  39  for  various  values  of  the  proportional  navigation  ratio.  From  the  same 
simulation,  it  is  possible  to  extract  similar  information  relating  to  the  miss  distance  due  to  a 
heading  error  (output  2  in  the  figure)  and  the  miss  distance  due  to  a  step  in  target  position 
(output  1  in  Figure  38).  These  are  presented  in  Figures  40  and  41  below.  It  is  clear  from  these 
figures  that  the  miss  distance  is  generally  smaller  as  the  navigation  ratio  increases. 
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Figure  39:  Miss  due  to  Target  Manoeuvre  for  Various  N' 


Flight  Time  (s) 


Figure  40:  Miss  due  to  Heading  Error  for  Various  N' 
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Figure  41:  Miss  due  to  Target  Displacement  for  Various  N' 


7.  Higher  Order  Guidance  System  Dynamics 


It  has  been  shown  that  when  the  guidance  system  dynamics  may  be  characterised  by  a  single 
time  constant,  then  it  is  possible  to  derive  closed  form  solutions  for  the  miss  distance. 
However,  it  is  known  from  reference  [1]  that  the  single  time  constant  representation  of  the 
guidance  system  seriously  underestimates  the  miss  distance.  For  a  more  realistic 
representation,  reference  [1]  recommends  the  use  of  a  fifth  order  transfer  function 
representation.  In  this  case,  no  closed  form  solution  is  possible.  Thus  the  only  option  is  to  use 
numerical  simulation. 


In  this  section,  the  fifth  order  binomial  model,  as  recommended  by  Zarchan  [1],  has  been 
adopted.  The  model  may  be  expressed  in  the  form. 


nL  =  N'Vcs 
A  (l  +  sr/5)5 


(88) 


Using  the  above  model,  it  is  shown  that  both  linear  models,  namely.  Model  A  and  Model  B, 
still  produce  equivalent  simulation  results  provided  that  the  correction  term  is  applied  to 
Model  A  as  previously  proposed  for  the  single  time  constant  case. 
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Firstly,  Model  A  (see  Figure  23)  is  modified  to  reflect  the  implementation  of  the  fifth  order 
guidance  system  representation  expressed  in  equation  (88)  above.  After  a  block  diagram 
restructure,  this  new  model  is  designated  as  Model  A5  and  is  shown  in  Figure  42  below.  Note 
that,  as  proposed  earlier  in  the  single  time  constant  case,  an  initial  condition  x  on  the  first 

integrator  downstream  of  the  entry  point  has  been  included  in  the  model  to  reflect  the  case  of 
a  step  in  target  displacement  y .  .  As  previously,  the  value  of  this  initial  condition  is  given  by 

the  expression. 


x. 


1C 


=  JL 


(89) 


MD  =  y(tF) 


Figure  42:  Block  Diagram  of  Model  A5 

Now,  following  the  adjoint  construction  rules  given  in  [1, 3],  it  is  an  easy  task  to  construct  the 
adjoint  of  Model  A5.  This  is  shown  in  Figure  43. 

In  a  similar  fashion,  the  alternative  adjoint  model  based  on  Model  B  is  constructed.  For  the 
fifth  order  system.  Model  B  is  designated  as  Model  B5.  This  is  shown  in  Figure  44  and  the 
corresponding  adjoint  model  is  shown  in  Figure  45. 

Both  models  were  then  programmed  in  Simulink  for  subsequent  adjoint  simulation.  The 
results  are  presented  in  Figure  46  through  to  Figure  48  inclusively.  It  is  clear  that  both  models 
yield  the  same  results  in  this  case  and  hence  any  of  the  models  may  be  used  interchangeably 
for  the  guidance  loop  analysis  provided  that  the  proposed  correction  strategy  derived  here  is 
adopted  when  using  Model  A. 
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Figure  43:  Adjoin  t  Block  Diagram  of  Model  A5 
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Figure  44:  Block  Diagram  of  Model  B5 
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Figure  45:  Adjoint  Block  Diagram  of  Model  B5 


Flight  Time  (s) 


Figure  46:  Miss  Distance  due  to  Target  Manoeuvre  for  5th  Order  System 
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Flight  Time  (s) 


Figure  47:  Miss  Distance  due  to  Heading  Error  for  5th  Order  System 


Figure  48:  Miss  Distance  due  to  Target  Displacement  for  5th  Order  System 
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8.  Conclusions 


In  this  report,  attention  is  focussed  on  the  models  and  tools  often  adopted  for  the  linear 
analysis  of  missile-target  engagements  in  two  dimensions.  Following  a  review  of  the 
non-linear  dynamic  equations  describing  the  problem,  two  linear  models  were  developed  for 
the  linear  analysis.  These  models  are  often  used  interchangeably  in  the  literature  for  missile 
guidance  loop  analysis.  However,  it  is  shown  here  that  under  certain  initial  conditions  of  the 
state  variables,  the  two  models  yield  significantly  different  results.  Hence,  caution  needs  to  be 
exercised  by  the  analyst  when  adopting  these  models  for  missile  homing  loop  analysis. 

Following  closer  analysis,  a  correction  factor  was  derived  which  may  be  used  to  render  the 
two  models  equivalent.  This  correction  factor  was  then  used  as  the  basis  for  deriving  closed 
form  solutions  to  the  problem  in  the  case  when  the  guidance  system  may  be  approximately 
represented  by  a  single  time  constant  system.  After  carrying  out  the  mathematics,  it  was 
demonstrated  that  both  models  yield  the  same  closed  form  solution  for  the  particular  inputs 
considered. 

The  models  were  then  implemented  in  Simulink  and  it  was  shown  how  they  could  be  used  for 
assessment  of  other  problems  in  which  the  closed  form  solutions  are  not  available.  In  addition 
to  insights  gained  on  system  performance,  another  benefit  of  having  derived  the  closed  form 
solutions  to  this  problem  was  the  observation  that  two  of  the  analytical  relationships  were 
indeed  connected.  Consequently,  it  was  shown  how  this  fact  could  be  utilised  for  the 
generation  of  performance  data  associated  with  the  study  of  multiple  targets  using  the  system 
response  to  a  heading  error  input. 

Finally,  the  power  and  flexibility  of  Simulink  was  utilised  to  demonstrate  how  easily  these 
linearised  engagement  models  and  their  adjoints  could  be  extended  to  accommodate  higher 
order  and  more  complex  dynamic  system  representations  of  the  guidance  loop. 


44 


UNCLASSIFIED 


UNCLASSIFIED 

9.  References 


DSTO-TR-2845 


1.  Zarchan,  P.,  "Tactical  and  Strategic  Missile  Guidance,"  5th  Ed.  Progress  in 
Astronautics  and  Aeronautics,  Vol.  124,  AIAA  Inc.  Washington  DC,  2007. 

2.  Shneydor,  N.A.,  "Missile  Guidance  and  Pursuit,  Kinematics,  Dynamics  and 
Control,"  Horwood  Publishing  Ltd,  Chichester,  UK,  1998. 

3.  Bucco,  D.,  "Aerospace  Applications  of  Adjoint  Theory,"  Monograph,  DSTO 
Fellowship  Program,  Australia,  Jan  2010. 

4.  Gutman,  S.,  "Applied  Min-Max  Approach  to  Missile  Guidance  and  Control," 
Progress  in  Astronautics  and  Aeronautics,  Vol.  209,  AIAA  Inc.  Washington  DC,  2005. 

5.  Weiss,  M.,  "Adjoint  Method  for  Missile  Performance  Analysis  on  State-Space 
Models,"  Journal  of  Guidance,  Control  and  Dynamics,  Vol.  28,  No.  2,  pp.  236-248, 
2005. 

6.  Taylor,  J.H.,  Price,  C.F.,  Siegel,  J.  &  Gelb,  A.,  "Covariance  Analysis  of  Nonlinear 
Stochastic  Systems  via  Statistical  Linearization,"  ASME,  Nonlinear  System  and 
Synthesis,  Vol.  2,  Techniques  and  Applications,  Feb  1981. 


UNCLASSIFIED 


45 


UNCLASSIFIED 

DSTO-TR-2845 

Appendix  A:  Non-linear  Engagement  Model  in 

MATLAB/Simulink 

A.l.  Details  of  the  Simulink  Model 

A.1.1  The  Simulink  Subsystem  Models 

The  missile  dynamics  equations  and  the  target  dynamic  equations  are  implemented  in 
Simulink  as  subsystems,  as  shown  in  the  figures  below. 


Figure  Al.  Missile  Dynamics  Subsystem 

The  Missile  Dynamics  Subsystem  (Figure  Al)  contains  all  of  the  blocks  required  to 
simulate  the  missile  dynamics  as  represented  by  the  non-linear  equations  given  in  the  text. 

In  a  similar  fashion,  the  Target  Dynamics  Subsystem  (Figure  A2)  contains  all  of  the  blocks 
required  to  simulate  the  target  dynamics  in  accordance  with  the  target  equations  given  in 
the  text.  In  order  to  accommodate  a  sudden  change  in  target  position,  a  switch  block  has 
been  included  in  the  Simulink  model.  The  switch,  in  this  case,  is  used  to  simulate  the 
different  values  of  the  target  position  to  be  represented  in  the  simulation  at  the  moment  of 
seeker  resolution  in  a  multiple  target  scenario.  It  was  necessary  to  insert  a  unit  delay  in  the 
loop  in  order  to  avoid  an  inherent  algebraic  loop  condition  causing  the  simulation  to  fail. 
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Switch 


Figure  A2.  Target  Dynamics 

The  implementation  of  the  relative  dynamics  is  shown  in  Figure  A3. 


Figure  A3.  Relative  Dynamics  Subsystem 
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Figure  A4.  DelR  Subsystem 


Figure  A5.  DelV  Subsystem 
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Figure  A6.  Closing  Vel  Subsystem 


Figure  A7.  Stop  Condition 


Figure  A8.  Lamda  Dot  Subsystem 
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Note  that  the  Relative  Dynamics  subsystem  block  shown  in  Figure  A3  is  composed  of  the 
subsystems  given  in  Figures  A4,  A5,  A6  and  A8.  The  simulation  stopping  condition  model  is 
shown  in  Figure  A 7. 


A.  1.2  The  Associated  MATLAB  Code 

The  missile-target  engagement  model  represented  in  Simulink  is  controlled  by  a  Graphical 
User  Interface  underpinned  by  MATLAB  code.  The  complete  MATLAB  code  is  listed  over  the 
following  pages. 


C  Editor  -  C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 


File 

Edit 

Text  Go  Cell  Tools  Debug  Desktop  Window  Help 

a 

a 

1 

F]  function  TwoD_EngageGUI 

“al 

2 

3  %  GUI  for  2D  Engagement  Simulation 

3 

% - 

4 

c 

-  %  Initialisation 

3 

6 

%  Get  screen  size  and  determine  figure  position  to  centre  figure 

7 

%  on  screen 

8 

- 

screenSize  =  get(Of  ' ScreenSize ' ) ; 

9 

- 

figureSize= [600  350]; 

10 

- 

figurePos  =  [ 

11 

(screenSize (3)  -  f igureSize (1) ) /2 

12 

(screenSize (4)  -  figureSize (2) ) /2 

13 

figureSize (1) 

14 

figureSize (2) ] ; 

15 

% 

16 

%  Default  parameters  for  Target  model 

17 

% 

18 

%  Target  Position  (m) 

19 

- 

appData.XTO  =  6500; 

20 

- 

appData.YTO  =  350; 

— 

21 

%  Target  flight  path  angle 

ZI 

22 

- 

appData.BetaO  =  0; 

zz 

23 

%  Target  Manoeuvre  (g) 

= 

24 

- 

appData.NT  =3; 

25 

%  Target  Speed  (m/s) 

26 

- 

appData.VT  =  300; 

27 

% 

28 

%  Default  parameters  for  Missile  Model 

29 

% 

TwoD EngageGUI  [In 

Col  20 

Figure  A9.  MATLAB  Code 
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6  Editor  -  C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 


[ext  Go  Cell  Tools  Debug  Desktop  Window  Help 

X 

! 

30 

%  Missile  Position  (ra) 

31  - 

appData.XMO  =  0; 

32  - 

appData. YMO  =  350; 

33 

%  Missile  Speed  (m/s) 

34  - 

appData. VM  =  1000; 

35 

%  Missile  Time  Constant 

36  - 

appData. Tau  =  .3; 

37 

%  Effective  Navigation  Ratio 

38  - 

appData.Np  =3; 

39 

%  Heading  Error  (degrees) 

40  - 

appData. HE  =0; 

41 

%  Time  to  go  when  target  resolution  occurs 

42  - 

appData. Thom  =1; 

43 

%  Target  Displacement  when  resolution  occurs 

44  - 

appData. Displace  =0; 

45 

% 

46 

%  Initial  Miss  Value 

47  - 

appData.Tmiss  =0; 

48  - 

appData. Miss  =0; 

49 

%  -  MAIN  FIGURE  - 

50  - 

GUIfig  =  f igure (  'Tag',  'GUIfig',  ... 

z: 

51 

■Units',  'pixels',  ... 

= 

52 

■Position',  figurePos,  ... 

= 

53 

■Name',  '2D  Engagement',  ... 

54 

'MenuBar',  'none',  ... 

55 

' Resize ' , ' off ' , . . . 

56 

' NumberTitle ' ,  'off',  ... 

57 

'Color',  get(0, ' Def aultUicontrolBackgroundColor ' ) ) ; 

58 

%  -  SUB  PANEL  -  Target  - 

V 

TwoD EngageGUI  Ln  34  Col  19 


Figure  A10.  MATLAB  Code  (continued) 


Figure  All.  MATLAB  Code  (continued) 
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B  Editor 

C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 

1  File  Edit 

[ext  Go  Cell  Tools  Debug  Desktop  Window  Help 

X 

O  1 

,  11  i  Ij,  •  ■  • 

□ 

88 

'String',  appData. YTO,  ... 

89 

•  Callback: • ,  @YTval_Callback)  ; 

90  - 

uicontrol (  ' Parent* ,  TgtPanel,  ... 

91 

’ Tag • ,  ’ YTvalTXT ' ,  ... 

92 

'Style' ,  'text' ,  ... 

93 

'Units',  'pixels',  ... 

94 

'Position',  [180  80  30  20],  ... 

95 

'String',  ' YT0 ’ ) ; 

96 

%  BetaO  value 

97  - 

uicontrol (  'Parent',  TgtPanel,  ... 

98 

'Tag' ,  'BetaO' ,  ... 

99 

'Style' ,  'edit' ,  ... 

100 

'Units',  'pixels',  ... 

101 

•Position',  [100  40  40  20] ,  ... 

102 

’ BackgroundColor ' ,  [111],  ... 

103 

'String',  appData . BetaO,  ... 

104 

'Callback',  @Beta0  Callback); 

105  - 

uicontrol (  'Parent',  TgtPanel,  ... 

= 

106 

'Tag',  'BetaOTxt',  ... 

II 

107 

’ Style ' ,  ' text ' ,  ... 

= 

108 

'Units',  'pixels',  ... 

109 

'Position',  [60  35  30  20],  ... 

110 

'String',  'BetaO'); 

= 

111 

%  Target  Speed 

112  - 

uicontrol (  'Parent',  TgtPanel,  ... 

113 

'Tag' ,  'VTval' ,  ... 

114 

'Style',  'edit',  ... 

115 

'Units',  'pixels',  ... 

116 

•Position',  [100  10  40  20] ,  ... 

v 

TwoD EngageGUI  Ln  96 

Col  16 

Figure  All.  MATLAB  Code  (continued) 


E  Editor  -  C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 


^il^Tdit 

Text  Go  Cell  Tools  Debug  Desktop  Window  Help 

a 

a 

118 

'String',  appData. VT,  ... 

A 

119 

■  Callback ' ,  @VTval_Callback) ; 

120  - 

uicontrol (  'Parent',  TgtPanel,  ... 

121 

'Tag',  ’ VTvalTxt ' ,  ... 

122 

' Style ' ,  ' text ' ,  ... 

123 

'Units',  'pixels',  ... 

124 

'Position',  [60  8  20  20],  ... 

125 

' String ' ,  ' VT ' ) ; 

126 

%  -  SUB  PANEL  -  Missile  - 

127  - 

MslPanel  =  uipanel (  'Parent',  GUIfig,  ... 

128 

'Tag',  'MslPanel',  ... 

129 

'Units',  'pixel',  ... 

130 

•Position',  [320  180  250  150],  ... 

131 

'Title',  'Missile  Parameters’); 

132 

%  -  Objects  in  Missile  Panel  - 

133 

%  XM0  Value 

134  - 

uicontrol (  'Parent',  MslPanel,  ... 

135 

' Tag ' ,  ’ XMval ’ ,  ... 

- 

136 

' Style ' ,  ' edit ' ,  ... 

= 

137 

'Units',  'pixels',  ... 

= 

138 

'Position',  [20  100  70  20] ,  ... 

— 

139 

'BackgroundColor',  [111],  ... 

— 

140 

'String',  appData.XMO,  ... 

= 

141 

•  Callback' ,  @XMval_Callback) ; 

— 

142  - 

uicontrol (  'Parent',  MslPanel,  ... 

143 

■Tag',  ’ XMvalTXT ' ,  ... 

144 

' Style ' ,  ' text ' ,  ... 

- 

145 

■Units',  'pixels',  ... 

146 

'Position',  [40  80  30  20],  ... 

1  4" 

1  .Sr.ri  nr* '  .  'XMOM; 

v 

_ 

TwoD EngageGUI  Ln  96 
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Figure  A13.  MATLAB  Code  (continued) 
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C  Editor 

C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 

□0® 

File  Edit 

[ext  Go  Cell  Tools  Debug  Desktop  Window  Help 

X 

□ 

148 

%  YMO  Value 

A 

149  - 

uicontrol!  'Parent',  MslPanel,  ... 

150 

'Tag' ,  ' YMval' ,  ... 

151 

'Style',  'edit',  ... 

152 

'Units',  'pixels’,  ... 

153 

'Position',  [160  100  70  20],  ... 

154 

■ 3ackgroundColor ' ,  [111],  ... 

155 

'String',  appData.YMO,  ... 

156 

•Callback:',  @YMval  Callback); 

157  - 

uicontrol (  'Parent',  MslPanel,  ... 

ise 

'Tag',  ' YMvalTXT ’ ,  ... 

159 

'Style' ,  'text' ,  ... 

160 

'Units',  'pixels',  ... 

161 

'Position',  [180  80  30  20],  ... 

162 

'String',  'YMO'); 

163 

%  Nprime  value 

164  - 

uicontrol(  'Parent',  MslPanel,  ... 

165 

'Tag',  'Nprime',  ... 

= 

166 

'Style' ,  'edit' ,  ... 

= 

167 

'Units',  'pixels',  ... 

= 

168 

■Position',  [60  40  40  20] ,  ... 

“ 

169 

' 3ackgroundColor ' ,  [111],  ... 

170 

'String',  appData.Np,  ... 

171 

'Callback',  @Nprime  Callback); 

172  - 

uicontrol!  'Parent',  MslPanel,  ... 

173 

'Tag',  'NprimeTxt',  ... 

174 

' Style ' ,  ' text ' ,  ... 

— 

175 

'Units',  'pixels',  ... 

176 

'Position',  [20  35  20  20] ,  ... 

V 

TwoD EngageGUI  [Tn  157 
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Figure  Aid.  MATLAB  Code  (continued) 
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C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 

□0® 

File  Edit 

Text  Go  Cell  Tools  Debug  Desktop  Window  Help 

-x  ? 

X 

176 

■Position',  3b  ZU  ZUJ,  ... 

, 

177 

' String* ,  *N' ' ' ) ; 

A 

178 

%  Missile  Speed 

179  - 

uicontrol!  'Parent',  MslPanel,  ... 

180 

'Tag',  'VMval',  ... 

181 

'Style' ,  'edit' ,  ... 

182 

'Units',  'pixels',  ... 

183 

'Position',  [60  10  40  20] ,  ... 

184 

'3ackgroundColor ' ,  [111],  ... 

185 

'String',  appData.VM,  ... 

186 

' Callback' ,  @VMval_Callback) ; 

187  - 

uicontrol!  'Parent',  MslPanel,  ... 

188 

'Tag',  'VMvalTxt',  ... 

189 

'Style' ,  'text' ,  ... 

190 

'Units',  'pixels',  ... 

191 

'Position',  [20  8  20  20],  ... 

192 

'String',  'VM'); 

193 

%  Guidance  Loop  Time  Constant 

— 

194  - 

uicontrol!  'Parent',  MslPanel,  ... 

= 

195 

'Tag' ,  'Tauval' ,  ... 

= 

196 

'Style' ,  'edit' ,  ... 

= 

197 

'Units',  'pixels',  ... 

— 

198 

•Position',  [180  40  40  20] ,  ... 

199 

■BackgroundColor ' ,  [111],  ... 

200 

'String',  appData.Tau,  ... 

201 

' Callback' ,  @GLTC_Callback) ; 

202  - 

uicontrol!  'Parent',  MslPanel,  ... 

203 

' Tag ' ,  ' TauvalTxt ' ,  ... 

204 

'Style' ,  'text' ,  ... 

205 

'Units',  'pixels',  ... 

V 

TwoD EngageGUI  Ln  186 
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Figure  A15.  MATLAB  Code  (continued) 
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£  Editor 

C:\Work  3\Reports\Two  Target\SimuIations\TwoD_EngageGUI.m 

0CS® 

File  Edit 

[ext  Go  Cell  Tools  Debug  Desktop  Window  Help 

•51  ? 

X 

206 

'Position',  [140  35  20  20],  ... 

A 

. 

207 

'String',  'Tau'); 

208 

%  -  SUB  PANEL  -  Sim  Inputs  - 

209  - 

SimPanel  =  uipanel (  'Parent',  GUIfig,  ... 

210 

'Tag',  'SimPanel',  ... 

211 

'Units',  'pixel',  ... 

212 

■Position',  [30  20  250  150],  ... 

213 

' ForegroundColor ' ,  [0.251  0  0.251],  ... 

214 

'Title',  'Simulation  Inputs’); 

215 

216 

217 

%  Target  Manoeuvre 

218  “ 

uicontrol (  'Parent',  SimPanel,  ... 

219 

'Tag',  'NTval ' ,  ... 

220 

'Style',  'edit',  ... 

221 

'Units’,  'pixels',  ... 

222 

■Position',  [120  104  40  20] ,  ... 

223 

'BackgroundColor ' ,  [111],  ... 

“ 

224 

'String',  appData.NT,  ... 

= 

225 

'  Callback • ,  @NTval_CallbacJc)  ; 

= 

226  - 

uicontrol (  'Parent',  SimPanel,  ... 

227 

'Tag',  'NTvalTxt',  ... 

228 

'Style* ,  'text' ,  ... 

= 

229 

'Units',  'pixels',  ... 

~~ 

230 

'Position',  [20  100  90  20] ,  ... 

231 

'String',  'NT  (g''s)'); 

232 

%  Heading  Error 

— 

233  - 

uicontrol (  'Parent',  SimPanel,  ... 

234 

'Tag',  'HEval',  ... 

-  :  : 

ssl 

1 

|  TwoD.EngageGUI  |Ln  214 
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Figure  A16.  MATLAB  Code  (continued) 


£  Editor 

C:\Work  3\Reports\Two  Target\SimuIations\TwoD_EngageGUI.m 

00® 

1  File  Edit 

Text  Go  Cell  Tools  Debug  Desktop  Window  Help 

-51  ? 

X 

-  t 

iay ■ ,  ■ ntvdi  ,  ... 

□ 

235 

'Style' ,  'edit' ,  ... 

A 

236 

'Units',  'pixels',  ... 

237 

•Position',  [120  73  40  20]  ,  ... 

238 

'BackgroundColor’,  [111],  ... 

239 

'String',  appData.HE,  ... 

240 

' Callback ' ,  @HEval_Callback) ; 

241  - 

uicontrol  (  'Parent',  SimPanel,  .  .  .| 

242 

'Tag',  'HEvalTxt',  ... 

243 

'Style' ,  'text' ,  ... 

244 

'Units',  'pixels',  ... 

245 

•Position',  [20  70  90  20] ,  ... 

246 

'String',  'HE  (deg)’); 

247 

%  Target  Displacement 

248  - 

uicontrol  (  'Parent',  SimPanel,  ... 

249 

■Tag' ,  ’ Disval ’ ,  ... 

250 

'Style' ,  'edit' ,  ... 

251 

'Units',  'pixels',  ... 

_ 

252 

'Position',  [120  42  40  20] ,  ... 

= 

253 

'BackgroundColor ' ,  [111],  ... 

254 

'String',  appData . Displace,  ... 

_ 

255 

■Callback' ,  @Disval_Callback) ; 

= 

256  - 

uicontrol (  'Parent',  SimPanel,  ... 

= 

257 

'Tag',  'DisvalTxt',  ... 

= 

258 

' Style ’ ,  ' text ' ,  ... 

259 

'Units',  'pixels',  ... 

260 

•Position',  [20  40  90  20] ,  ... 

261 

'String',  'Displace  (m)  ' ) ; 

262 

%  Time  at  which  target  resolution  occurs 

263  - 

uicontrol  (  'Parent',  SimPanel,  ... 

,V| 

TwoD EngageGUI  Ln  241 
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Figure  A17.  MATLAB  Code  (continued) 
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Figure  A18.  MATLAB  Code  (continued) 


C  Editor 

C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 

1  File  Edit 

Text  Go  Cell  Tools  Debug  Desktop  Window  Help 

-si  l* 

X 

r  ,  [jju  ju  iuj  ,  .  .  . 

■■ 

3 

292 

'BaclcgroundColor ' ,  [111],  ... 

* 

293 

'String',  appData.Tmiss,  ... 

294 

'Callback',  @Tmissval_Callback) ; 

295  - 

uicontrol (  'Parent',  GUIfig,  ... 

296 

'Tag',  'TmissTxt',  ... 

297 

•Style' ,  'text' ,  ... 

298 

'Units',  'pixels’,  ... 

299 

•Position',  [348  95  50  20] ,  ... 

300 

'String',  'Time  ( s ) ’ ) ; 

301 

%  Miss  Distance 

302  - 

uicontrol  (  'Parent',  GUIfig,  .  .  .| 

303 

'Tag',  'Missval',  ... 

304 

'Style' ,  'edit' ,  ... 

305 

'Units',  'pixels',  ... 

306 

•Position’,  [450  120  60  20] ,  ... 

307 

'BackgroundColor ' ,  [111],  ... 

308 

'String',  appData .Miss,  ... 

309 

'Callback' ,  @Missval_Callback) ; 

= 

310  - 

uicontrol (  'Parent',  GUIfig,  ... 

ii 

311 

'Tag',  'MissTxt',  ... 

= 

312 

'Style' ,  'text' ,  ... 

Z 

313 

'Units',  'pixels',  ... 

= 

314 

•Position',  [447  95  70  20] ,  ... 

= 

315 

'String',  'Miss  Dist  (m)  ' ) ; 

316 

317 

% - End  of  Parameter  Entry  For  GUI - 

318 

%  Save  application  data 

319  - 

guidata (GUIfig, appData) ; 

320  - 

L  end 

V 

TwoD EngageGUI  Ln  302 
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Figure  A19.  MATLAB  Code  (continued) 
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Figure  A20.  MATLAB  Code  (continued) 


Figure  All.  MATLAB  Code  (continued) 
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Figure  All.  MATLAB  Code  (continued) 


C  Editor 

C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 

1  File  Edit 

Text  Go  Cell  Tools  Debug  Desktop  Window  Help 

X 

407  - 

appData  =  guidata (hObject) ; 

a| 

408  - 

handles  =  guihandles (hObject) ; 

409  - 

HdgErr  =  str2num(get (handles .HEval, ' String' ) ) ; 

410  - 

appData. HE  =  HdgErr; 

411  - 

3et (handles .HEval, ' string ' , num2str (HdgErr) ) ; 

412  - 

guidata (hObject, appData) ; 

413  - 

^  end 

414 

%  Target  Displacement  Callback 

415 

E 

! function  DisvalCallback (hObject, eventdata) 

416  - 

appData  =  guidata (hObject) ; 

417  - 

handles  =  guihandles (hObject) ; 

418  - 

Disp  =  str2num (get (handles. Disval, 'String' )) ; 

419  - 

appData. Displace  =  Disp; 

420  - 

set (handles .Disval, ' string* , num2str (Disp) ) ; 

421  - 

guidata (hObject, appData) ; 

422  - 

end] 

423 

%  Homing  Time  Callback 

424 

E 

function  Thomval  Callback (hObject, eventdata) 

= 

425  - 

appData  =  guidata (hObject) ; 

= 

426  - 

handles  =  guihandles (hObject) ; 

427  - 

Thom  =  strgnum (get (handles .Thomval, ’ String' )) ; 

428  - 

appData. Thom  =  Thom; 

= 

429  - 

set (handles .Thomval, ' string’ , num2str (Thom) ) ; 

= 

430  - 

guidata (hObject, appData) ; 

431  - 

L  end 

432 

%  Guidance  Loop  Time  Constant  Callback 

433 

E 

function  GLTC  Callback (hObject, eventdata) 

434  - 

appData  =  guidata (hObject) ; 

435  - 

handles  =  guihandles (hObject) ; 

v 

TwoD EngageGUI  /  Disval Callback  Ln  422 
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Figure  A23.  MATLAB  Code  (continued) 
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C  Editor  -  C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m  f»~|fD]fx] 


File 

Edit 

[ext  Go  Cell  Tools  Debug  Desktop  Window  Help 

BE 

□ 

436 

- 

Tau  =  str2num(get (handles .Tauval, ' String' )) ; 

A 

□ 

437 

- 

appData.Tau  =  Tau; 

438 

- 

set (handles .Tauval, ' string' , num2str (Tau) ) ; 

439 

- 

guidata (hObject, appData) ; 

440 

- 

end 

441 

442 

443 

% - 

444 

%  Simulation  Callback 

445 

% - 

446 

□  function  RunSim  Callback  (hObject,  eventdata) 

447 

- 

appData  =  guidata (hObject) ;  %% 

448 

- 

handles  =  guihandles (hObject) ;  %% 

449 

% 

450 

%  Script  to  run  the  2D  missile/target  engagement  using  Prop  Nav  guidance 

451 

% 

452 

%  Parameters 

453 

- 

Np  =  appData. Np;  %  Navigation  constant 

— 

454 

- 

HE  =  appData. HE/57. 3;  %  Heading  Error 

455 

- 

Thom  =  appData. Thom;  %  Time  to  go  for  step  in  tgt  displacement 

= 

456 

- 

Displace  =  appData. Displace;  %  Step  in  tgt  displacement 

— 

457 

- 

Tau  =  appData.Tau;  %  Msl  Guidance  Loop  time  constant 

~ 

458 

— 

NT  =  appData. NT* 9. 8;  %  Tgt  maneuver 

= 

459 

%  Initial  Target  Data 

— 

460 

- 

XTO  =  appData.XTO; 

■ 

461 

- 

YTO  =  appData.YTO; 

462 

- 

YT1  =  YTO  +  Displace; 

— 

463 

- 

VT  =  appData.VT; 

464 

- 

BetaO  =  appData. BetaO/57. 3; 

465 

%  Initial  Missile  Data 

|vj| 

TwoD EngageGUI  /  Disval Callback  Ln  422 
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Figure  AM.  MATLAB  Code  (continued) 


Q  Editor 

C:\Work  3\Reports\Two  Target\Simulations\TwoD_EngageGUI.m 

1  File  Edit 

lext  Go  Cell  Tools  Debug  Desktop  Window  Help 

X 

466  - 

XMO  =  appData.XMO; 

□ 

467  - 

YMO  =  appData. YMO; 

468  - 

VM  =  appData.VM; 

469 

%  Relative  Geometry 

470  - 

DelXO  =  XTO  -  XMO; 

471  - 

DelYO  =  YTO  -  YMO; 

472  - 

LamdaO  =  atan2 (DelYO, DelXO) ; 

473  - 

L  =  asin  (VT*sin (BetaO+LamdaO) /VM) ; 

474  - 

ThetaO  =  LamdaO  +  L; 

475  - 

VMXO  =  VM*cos (ThetaO  +  HE); 

476  - 

VMYO  =  VM*sin (ThetaO  +  HE); 

477 

%  Load  up  these  initial  parameters  into  the  Simulink  model 

478 

%  make  sure  the  Simulink  model  is  open 

479  - 

SimMod  =  'TwoDMod'; 

480 

%  set  time  constant) 

481  - 

Den  =  [’['  num2str(Tau)  '  1]']; 

482  - 

set_param( [SimMod  ' /AP ' ] , 'Numerator ' , ' [1] ' ,  . . . 

483 

' Denominator ’ ,  Den) ; 

= 

484 

%  set  tgt  man  value 

= 

485  - 

set  param( [SimMod  ' /NT ' ] , ' Value ' , num2str (NT) ) ; 

486 

%  set  parameters  in  Target  Dynamics  subsystem 

” 

487  - 

set  param( [SimMod  '/Target  Dynamics/X' ] , ’ InitialCondition' , num2str (XTO) ) ; 

488  - 

set  param( [SimMod  '/Target  Dynamics/Beta'],... 

= 

489 

' InitialCondition' , num2str (BetaO) ) ; 

490  - 

set_param( [SimMod  '/Target  Dynamics/YTO 'Value ', num2str (YTO) ) ; 

491  - 

set_param( [SimMod  '/Target  Dynamics/YTl '], 'Value ', num2str (YT1) ) ; 

492  - 

set_param( [SimMod  '/Target  Dynamics/Thom' ],' const ’, num2str (Thom) ) ; 

— 

493  - 

set  param( [SimMod  '/Target  Dynamics/VTO 'Value ’, num2str (VT) ) ; 

494 

%  set  parameters  in  Missile  Dynamics  subsystem 

v| 

...D EngageGUI / RunSim Callback  Ln  480 
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Figure  A25.  MATLAB  Code  (continued) 
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Figure  A26.  MATLAB  Code  (continued) 
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Appendix  B:  Closed  Form  Solutions  for  Single  Time 

Constant  Guidance  System 


If  the  guidance  system  can  be  approximated  by  a  single  time  constant  parameter  with  transfer 
function  of  the  form,  (see  equation  (25)  in  the  text) 

N'  (Bl) 

5(1  +  zs)  ’ 

then,  for  integral  values  of  TV  ,  it  has  been  shown  in  the  text  how  to  derive  useful  closed  form 
solutions  to  the  linear  engagement  problem.  In  particular,  typical  closed  form  solutions  have 
been  provided  therein  for  the  specific  case  when  N  =  3 .  These  analytical  expressions  are 
useful  when  producing  design  curves  during  the  preliminary  stages  of  a  guided  missile 
program.  In  this  appendix,  closed  form  solutions  for  integral  values  of  N  in  the  range  [4  -  6] 
are  provided. 


B.l.  Solution  for  N'  =  4 


MD 


n „ 


t 2  , , 

t 

=  ^eF 

1-^ 

2 

N'=4 

3r_ 

(132) 


mdha 

HE 


r  t  e  1 

=  -VMtFe"rlT 

1 

b1 

+ 

M  F 

N'=4 

1 

1 _ 

(B3) 


MDn 

YT 


[  t  e  1 

Ci 

II 

1 

+ 

if 

_  T  6r  _ 

(B4) 


B.2.  Solution  for  N'  =  5 


MD, 


n., 


1  t2  . 

r  2c  t2  1 

C 

11 

1 

+ 

•> 

2 

N'=5 

L  3r  12r  J 

(B5) 
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MD 


_ HA  | 

HE 


r  3 tF 

tl 

tl 

=  -VMtFe'rlz 

1  _  F 

+  F 

F 

M  F 

N'-5 

2  T 

2r2 

24r3  _ 

(B6) 


MD, 


YT 


=  e 


1  fF 

1 — -  +  F 


2  E 

lF 


2t  It2  24 H 


(B  7) 


B.3.  Solution  for  N'  =  6 


MD, 


n. 


-  —  e'r'T 


t  tl  ti 

F  .  F  F 


i  -m+ 
t  4r2 


60r3 


MD 


HA 

HE 


=  -V  tFe~ 

M  F 


2  A  tl  tl 


tl 


1 - ^  +  -J- - ^  + 

t  t2  6r3  120r4 


MDn 

YT 


=  e 


2 1 
T  T 


F  i  ~F 


tl  tl  tl 

+  ^ - ^  +  - - - 

2  6r3  120r4 


(B8) 


(B9) 


(BIO) 


The  above  closed  form  solutions  have  been  verified  against  the  Simulink  models. 

It  is  interesting  to  note  from  the  above  expressions  that  the  miss  distance  due  to  a  target 
displacement  appears  to  be  directly  related  to  the  miss  distance  due  to  a  heading  error  in  the 
following  way. 


MD „ 
YT 


N'-n 


1  MD, 


tF  (■ -VmHE ) 


(Bll) 


where  n  is  an  integer  greater  than  2.  In  fact,  through  simulation  of  the  adjoint  models  in 
Simulink,  it  was  straightforward  to  verify  that  this  relationship  holds  true  even  for  non¬ 
integer  values  of  N  and  also  for  higher  order  systems  such  as  the  fifth  order  binomial 
system.  Further  investigations  of  this  relationship  are  carried  out  in  Appendix  C. 
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Appendix  C:  Alternative  Approach  to  Performance 
Prediction  of  Miss  Distance  due  to  a  Step  in  Target 

Displacement 

The  relationship  given  in  equation  (1311)  expresses  the  fact  that  the  miss  distance  due  to  a  step 
in  target  displacement,  for  a  single  time  constant  system,  may  be  determined  by  knowing  the 
miss  distance  due  to  an  initial  error  in  heading  angle.  That  is,  imposing  an  initial  heading 
angle  on  the  interceptor  of  the  form 


HA  = 


YT 


(Cl) 


and  assuming  there  are  no  other  inputs  to  the  system,  then  the  miss  distance  output  for  the 
problem  is  that  due  to  a  step  in  target  displacement  of  magnitude  YT .  This  is  certainly  true  of 
the  closed  form  solutions  given  in  the  main  text  and  in  Appendix  B  above.  The  investigation 
here  will  focus  on  whether  this  is  true  for  non-integral  values  of  N  and  also  for  higher  order 
guidance  systems. 

Let  the  initial  heading  angle  condition,  equivalent  to  equation  (Cl)  above,  be  denoted  by 
HAyt.  Imposing  this  condition  on  the  first  integrator  in  Figure  Cl  is  tantamount  to  imposing  a 
step,  of  magnitude  YT,  in  target  displacement  condition  on  the  second  integrator  in  the 
diagram. 

C.l.  Linear  System  Investigations 

For  the  linear  analysis,  either  Model  A  or  Model  B  may  be  used  in  this  case  as  they  have  been 
shown  to  produce  equivalent  results  when  the  input  is  heading  error.  For  the  present 
investigation.  Model  A  is  employed.  The  Model  A  forward  loop  block  diagram  incorporating 
the  HAyt  initial  condition,  as  expressed  in  equation  (Cl),  is  presented  in  Figure  Cl. 


MD  =  y(tF) 


Figure  Cl:  Model  A  with  HAyt  condition 
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The  adjoint  model  of  Figure  Cl  may  be  easily  constructed  using  the  adjoint  construction  rules 
and  is  displayed  in  Figure  C2. 


Figure  C2:  Adjoint  of  Model  A  with  HAYt  condition 


The  adjoint  model  is  employed  here  because  of  the  desire  to  compare  the  miss  distance 
profiles  in  only  one  run.  The  Simulink  representation  of  the  adjoint  block  diagram  for  carrying 
out  the  adjoint  simulation  is  shown  in  Figure  C3  below. 


Minus  Np  Vc 


Figure  C3:  Simulink  Representation  of  Adjoint  Model  A 

In  this  case,  output  channels  1, 2  and  4  shown  in  Figure  C3  are  of  interest.  Output  channels  1 
and  4  provide  the  data  to  compute  the  miss  distance  response  to  a  step  in  target  displacement 
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of  magnitude  YT=yic  while  output  channel  2  provides  the  data  for  the  same  response  via  the 
indirect  approach  proposed  here  using  the  HAyt  condition. 

The  miss  distance  results  generated  for  the  case  of  a  single  time  constant  system  are  compared 
in  Figure  C4.  For  selected  values  of  the  navigation  ratio,  the  plots  compare  the  miss  distance 
due  to  a  target  displacement  (  YT=yic  )  derived  using  the  conventional  approach  with  that 
using  the  FIAyt  approach.  It  is  clear  from  the  figure  that  the  results  agree  even  if  the  value  of 
N  is  non-integral.  In  Figure  C5,  simulation  results  generated  on  the  basis  of  a  fifth  order 
guidance  system  are  presented.  Again,  the  FIAyt  approach  is  seen  to  agree  favourably  with 
the  conventional  approach  to  the  miss  distance  problem  for  this  case. 

In  conclusion,  for  a  linear  guidance  system,  the  miss  distance  performance  due  to  a  step  in 
target  displacement  may  be  determined  on  the  basis  of  an  appropriate  heading  error  input 
thus  eliminating  the  need  to  model  the  target  displacement  effect  conventionally.  This 
approach  also  obviates  the  need  to  incorporate  the  correction  factor  in  Model  A  as  outlined 
earlier  in  the  text. 


Figure  C4:  Comparison  ofHAYT  approach  with  direct  approach  for  single  time  constant  system 
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Flight  Time  (s) 


Figure  C5:  Comparison  of  HAyt  approach  with  direct  approach  for  fifth  order  system 


C.2.  Non-linear  System  Investigations 

In  this  section,  the  accuracy  of  the  HAyt  approach  is  investigated  in  the  context  of  the  non¬ 
linear  simulation  model.  Figure  C6  shows  the  typical  trajectories  involved  when 
implementing  either  the  conventional  approach  or  the  FIAyt  approach.  As  stated  earlier,  the 
conventional  approach  here  means  applying  the  target  displacement  initial  condition  at  a 
certain  time  to  go  before  intercept  while  the  HAyt  approach  means  applying  an  equivalent 
initial  condition  in  accordance  with  equation  (Cl).  In  order  to  implement  the  HAyt  condition, 
it  is  necessary  to  estimate  the  intercept  time  at  the  start  of  the  simulation.  In  this  regard,  the 
following  approximation  is  utilised,  namely,  tF  =  R/  V c  ■  In  Figure  C7  is  presented  the  miss 

distance  profile  generated  by  carrying  out  multiple  runs  of  the  non-linear  simulation  model 
employing  the  HAyt  initial  condition.  All  other  inputs  are  set  to  zero  in  this  case. 
Superimposed  on  this  plot  is  the  corresponding  miss  distance  profile  produced  by  the  linear 
model.  The  results  agree  favourably. 
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Figure  C6:  Comparison  of  Conventional  and  HAYt  based  Trajectories 
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Figure  C7:  Comparison  of  Linear  and  Non-linear  Results  using  the  HAyt  Condition 


Next,  the  comparison  of  the  HAyt  based  results  with  those  generated  using  the  conventional 
approach  is  made,  using  the  non-linear  model  outlined  in  Section  2  in  the  text.  All  inputs 
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except  the  initial  target  displacement  are  set  to  zero  in  this  case.  The  simulation  results  are 
shown  in  Figure  C8  below. 

The  non-linear  simulations  are  repeated  in  order  to  examine  the  effects  on  the  miss  distance 
results  due  to  other  input  conditions,  such  as  target  manoeuvre  and  heading  error.  The  total 
miss  distance  results  are  shown  in  Figure  C9  below.  Note  that  these  results  are  based  on  a 
target  displacement  of  60  m. 
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Figure  C8:  Non-linear  Results  Comparing  Conventional  and  HAyt  Approach  (Single  Input) 


Finally,  using  the  non-linear  simulation  model,  the  effect  of  different  target  displacements  on 
the  total  miss  distance  results  is  investigated.  The  plots  are  summarised  in  Figure  CIO  and 
compare  the  miss  distance  profiles  generated  using  the  conventional  approach  with  those 
generated  using  the  HAyt  approach.  It  is  clear  from  the  plots  that  the  proposed  HAyt 
approach  to  the  target  displacement  problem  leads  to  acceptable  results  in  this  case. 
Moreover,  for  flight  times  larger  than  1  second,  the  plots  indicate  that  the  PN  guidance  law 
works  well  in  eliminating  errors  due  to  the  initial  target  displacement  and  therefore  the  miss 
distance  is  insensitive  to  the  magnitude  of  the  initial  target  displacement.  However,  if  the 
homing  time  is  less  than  1  second,  the  miss  distance  only  reduces  as  the  magnitude  of  initial 
target  displacement  is  reduced. 


66 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2845 


0 

o 

c 

TO 

0 

b 


TO 

O 


Flight  Time  (s) 


Figure  C9:  Non-linear  Results  Comparing  Conventional  and  HAyt  Approach  (Multiple  Inputs) 
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Figure  CIO:  Comparison  of  HAyt  Method  with  Conventional  Method 
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